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If the primordial fluctuations are non-Gaussian, then this non-Gaussianity will be apparent in 
the cosmic microwave background (CMB) sky. With their sensitive all-sky observation, MAP and 
Planck satellites should be able to detect weak non-Gaussianity in the CMB sky. On large angular 
scale, there is a simple relationship between the CMB temperature and the primordial curvature 
perturbation: AT/T — — $/3. On smaller scales; however, the radiation transfer function becomes 
more complex. In this paper, we present the angular bispectrum of the primary CMB anisotropy that 
uses the full transfer function. We find that the bispectrum has a series of acoustic peaks that change 
a sign, and a period of acoustic oscillations is twice as long as that of the angular power spectrum. 
Using a single non-linear coupling parameter to characterize the amplitude of the bispectrum, we 
estimate the expected signal-to-noise ratio for COBE, MAP, and Planck experiments. In order 
to detect the primary CMB bispectrum by each experiment, we find that the coupling parameter 
should be larger than 600, 20, and 5 for COBE, MAP, and Planck experiments, respectively. Even 
for the ideal noise- free and infinitesimal thin-beam experiment, the parameter should be larger than 
, 3. We have included effects from the cosmic variance, detector noise, and foreground sources in 

the signal-to-noise estimation. Since the simple inflationary scenarios predict that the parameter 
' is an order of 0.01, the detection of the primary bispectrum by any kind of experiments should 

be problematic for those scenarios. We compare the sensitivity of the primary bispectrum to the 
primary skewness and conclude that when we can compute the predicted form of the bispectrum, it 
becomes a "matched filter" for detecting the non-Gaussianity in the data, and much more powerful 
tool than the skewness. For example, we need the coupling parameter of larger than 800, 80, 70, and 
' 60 for each relevant experiment in order to detect the primary skewness. We also show that MAP 

. and Planck can separate the primary bispectrum from various secondary bispectra on the basis of 

the shape difference. The primary CMB bispectrum is a test of the inflationary scenario, and also 
' a probe of the non-linear physics in the very early universe. 
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Why measure the bispectrum of the cosmic microwave background (CMB) radiation anisotropy? Simple inflationary 
models predict that the CMB anisotropy field is nearly random Gaussian, and that two-point statistics completely 
specify statistical properties of CMB. However, our universe may not be so simple. Higher order statistics, such as 
. , the three-point correlation function, or its harmonic transform, the angular bispectrum, are potential probes of the 
physics of generating the primordial fluctuations. Since gravitationally induced non-linearities are small at z ~ 1300, 
5-H \ CMB is expected to be the best probe of the primordial non-Gaussianity Q . 

In the inflationary scenario , the quantum fluctuations of the scalar (inflaton) field generate the observed matter 
and radiation fluctuations in the universe |^|-^| . In the stochastic inflationary scenario of Starobinsky ]l0|] , the quantum 
fluctuations decohere to generate the classical fluctuations. There are two potential sources of non-Gaussianity in this 
inflationary model: (a) the non-linear coupling between the classical inflaton field and the observed fluctuation field, 
and (b) the non-linear coupling between the quantum noise field and the classical fluctuation field. The former has 
been investigated by Salopek and Bond ]H| , while the latter has been explored by Gangui et al. [Q. Calzetta and 
Hu jl3| and Matacz [Q present an alternative treatment of the decoherence process that leads to different results 
for the primordial density perturbation from those obtained by Starobinsky [jlOfl - Matacz's treatment makes similar 
predictions for the level of non-Gaussianity to the Starobinsky's treatment [|14|| . These studies conclude that in the 
slow roll regime, the fluctuations are Gaussian. However, features in the inflaton potential can produce significant 
non-Gaussianity [fa. 
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There have been claims for both the non-detection [[l6| and the detection 0,^8) of the non-Gaussianity in the 
COBE map. Banday, Zaroubi and Gorski argued the non-cosmological origin of the COBE non-Gaussianity. 
MAP and Planck will measure the fluctuation field down to angular scales ~ 0?2 and 0?1, and test these claims. 

Previous work on the primary non-Gaussianity has focused on very large angular scale, where the temperature 
fluctuations trace the primordial fluctuations. This is valid on the COBE scale. For MAP and Planck; however, we 
need the full effect of the radiation transfer function. In this paper, we develop a formalism for doing this, and then 
present numerical results. Both the formalism and the numerical results are main results of this paper. We also 
discuss how well we can separate the primary bispectrum from various secondary bispectra. 

This paper is organized as follows. Sec. || defines the bispectrum, the Gaunt integral, and particularly the new 
quantity called the "reduced" bispectrum, which plays a fundamental role in estimating the physical property of the 
bispectrum. Sec. Ill formulates the primary bispectrum that uses the full radiation transfer function, and presents 
the numerical results of the primary bispectrum and the skewness. Sec. IV estimates the secondary bispectra from 



the coupling between the Sunyaev-Zel'dovich and the weak lensing effects [20|-|22j, and from the extragalactic radio 
and infrared sources. Sec. [v] studies how well we can measure each bispectrum, and how well we can discriminate 
among various bispectra. Sec. ^ is devoted to further discussion and our conclusion. 



II. DEFINING THE "REDUCED" BISPECTRUM 



The observed CMB temperature fluctuation field AT(n)/T is expanded into the spherical harmonics: 



a lm ee I d J &—^-Y^(&), (1) 



where hats denote unit vectors. The CMB angular bispectrum is given by 



hhh = \ a hrn 1 0.l 2 m 2 a l 3 m 3 ) , \4) 



and the angle-averaged bispectrum is defined by 
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B hhh ee V ( h h h ) B,"Y7 2m3 , (3) 



mi TO2 wis 



where the matrix is the Wigner-3j symbol. The bispectrum B™^ 2 " 13 must satisfy the triangle conditions and selection 
rules: mi + m,2 + m-3 = Q, h+h + h = even, and \k — lj \ < h <h + lj for all permutations of indices. Thus, B™l™ 2Tri3 
consists of the Gaunt integral, G^l™ 2 ™ 3 , defined by 



Gl2T m ^ J d 2 m imi (n)Y hm2 (h)Y hm3 (h) 



(2Ji + 1) (2l 2 + 1) (2Z 3 + 1) fix h h \ ( h h h 
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gmim 2 m 3 ■ j satisfies all the conditions mentioned above. 

~ J lll 2 L 3 ' 

Given the rotational invariance of the universe, Bi t i 2 i 3 is written as 

where bi 1 i 2 i 3 is an arbitrary real symmetric function of l\, and I3. This form of equation (^|) is necessary and 
sufficient to construct generic -BJ^J 12 ™ 3 under the rotational invariance. Thus, we shall frequently use &; 1 ; 2 ; 3 instead 
of B"Y Q " i2?Tl3 in this paper, and call this function the "reduced" bispectrum, as bi t i 2 i 3 contains all physical information 
in B^lJ 12 " 13 . Since the reduced bispectrum does not contain the Wigner-3j symbol that merely ensures the triangle 
conditions and selection rules, it is easier to calculate and useful to quantify the physical properties of the bispectrum. 
The observable quantity, the angle-averaged bispectrum Bi t i 2 i 3 , is obtained by substituting equation (|5|) into (|^), 



(2/x + l) (2i 2 + l) (2Z 3 + 1) ( h h h 
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where we have used the identity: 
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Alternatively, one can define the bispectrum in the flat-sky approximation, 

(o(li)o(li)a(l3)) - (2^) 2 <S (2) (l! + h + la) ^(li.la, 1 3 ), (8) 



where 1 is the two dimensional wave-vector on the sky. This definition of B(\\, 1 2 , 13) corresponds to equation (|^), 
given the correspondence of G™!™*™ 3 -> <5 (2) (h + 1 2 + I3) in the flat-sky limit Q. Thus, 



I2, 13) (flat-sky approximation), (9) 



is satisfied. This fact also would motivate us to use the reduced bispectrum bi t i 2 i 3 rather than the angular averaged 
bispectrum Bi x i 2 i 3 . Note that bi 1 i 2 i 3 is similar to i?; 1 z 2 ; 3 defined by Magueijo [Esl. The relation is &i 1 z 2 i 3 = \/HrBi x i 2 i 3 . 



III. PRIMARY BISPECTRUM AND SKEWNESS 
A. Model of the primordial non-Gaussianity 

If the primordial fluctuations are adiabatic scalar fluctuations, then 

__$(k) 5 „(fc)y4(k), (10) 

where $(k) is the primordial curvature perturbation in the Fourier space, and ffr; (fc) is the radiation transfer function. 
a; TO thus takes over the non-Gaussianity, if any, from $(k). Although equation (Wuj is valid only if the universe is flat, 
it is straightforward to extend this to an arbitrary geometry. The isocurvature fluctuations can be similarly calculated 
by using the entropy perturbation and the proper transfer function. 
In this paper, we explore the simplest weak non-linear coupling case: 

$(x) = 4> L (x) + f NL (*|(x) - ($| (x))) , (11) 

in real space, where denotes the linear gaussian part of the perturbation. ($(x)) = is guaranteed. Henceforth, 

we shall call Jnl the non-linear coupling constant. This model is based upon the slow-roll inflationary scenario. 
Salopek and Bond and Gangui et al. [Q found that /jvl is given by a certain combination of the slope and the 
curvature of the inflaton potential. In the notation of Gangui et al., $3 = 2/jvl- Gangui et al. found that <I>3 ~ 10 -2 
in the quadratic and the quartic inflaton potential models. 
In the Fourier space, $(k) is decomposed into two parts: 

$(k)=$ L (k)+$ A r i (k), (12) 

and accordingly, 

O-lrn = af m + a Wn, (13) 

where $jvi(k) is the non-linear part defined by 

d 3 p 



NL 



[ $ i (k + p)ci>2(p)-(2 7 r) 3 ^ 3 )(k)<$|(x)) 



(14) 



(2tt) 3 

One can confirm that ($(k)) = is satisfied. In this model, a non-vanishing component of the $(k)-field bispectrum 
is 

($ L (k 1 )$L(k 2 )$ JVL (k 3 )) = 2(2 7 r) 3 <5( 3 >(k 1 + k 2 + k 3 )/ A r L F $ (fc 1 )P <& (fc 2 ), (15) 

where P$(fc) is the linear power spectrum given by (^l0^i)^l(}^2)) = (27r) 3 P$(fci)(5^ 3 ^(ki + k 2 ). We have also used 
($ L (k + p)$*(p)) = (2n) 3 P*(p)S^(k), and ($|(x)) =Wr 3 / d 3 kP„(k). 

Substituting equation ( [To|) into (g), using equation ( [l5| ) for the $(k)-field bispectrum, and then integrating over 
angles kj., ka, and k3, we obtain the primary CMB angular bispectrum, 
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where 

6f(0 = - / tfdkP^gTiimikr), (17) 

JO 

6f L (r) = - / k 2 dkf NL 9Ti(k)ji{kr). (18) 

Note that b\ (r) is a dimensionless quantity, while bf L (r) has a dimension of L -3 . 

One confirms that the form of equation @ holds. Thus, the reduced bispectrum, b^i^ — S™^™ 2 ™ 3 (^™/ 2 " 
(Eq.(^|)), for the primordial non-Gaussianity is 



m 2 m 3 \ 
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y prirnary = 2 Odr [6^ (0&f 2 (06^ (0 + ^ (0^ W^s M + ^ W&f 2 (0^ 3 W] ■ (19) 

Jo 

tfil™3 TV 1S ^ uu y specified by a single constant parameter /atl, as the cosmological parameters will be precisely 
determined by measuring the CMB angular power spectrum Ci (e.g., |24|]). It should be stressed again that this is the 
special case in the slow-roll limit. If the slow-roll condition is not satisfied, then Jnl = fNL{ki, &2, ^3) a t equation (|l5| ) 
|pq] . Wang and Kamionkowski |25fl have developed the formula to compute B; 1 z 2 i 3 from the generic form of $(k)-ficld 
bispectrum. Our formula (Eq.Ql6|)) agrees with theirs, given our form of the <f>(k)-field bispectrum (Eq.([l5|)). 

Even if the inflation produced Gaussian fluctuations, Pyne and Carroll pointed out that the general relativistic 
second-order perturbation theory would produce terms of /jvl ~ O(l) P^| . For generic slow-roll models, these terms 
dominate the primary non-Gaussianity. 
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FIG. 1. This figure shows bf(r) (Eq.@) and bf L (r) (Eq.@), the two terms in our calculation of the primary CMB 
angular bispectrum, as a function of r. Various lines in the upper panel show + l)&; L (r)/27r] x 10 10 , where r = c (to — t), 
at t = 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, and 1.6 x t* (decoupling time), while [bf L (r)f~l] x 10 10 are shown in the lower panel, to is 
the present-day conformal time. Note that cto = 11.8 Gpc, and ct» = 235 Mpc in our cosmological model chosen here. The 
thickest solid line in the upper panel is the CMB angular power spectrum [1(1 + l)C;/27r] x 10 10 . Ci is shown for comparison. 
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B. Numerical results of the primary bispectrum 



We evaluate the primary CMB bispectrum (Eqs. (|T6|)— (]Tg|)) numerically. We compute the full radiation transfer 
function griik) with the CMBFAST [^7| code, and assume the single power law spectrum, P$(fc) oc fc™~ 4 , for the 
primordial curvature fluctuations. The integration over k (Eqs.(|l7|) and ([i"8|)) is done by the algorithm used in 
CMBFAST. The cosmological model is the scale-invariant standard cold dark matter model with fl m = 1, = 0, 
fib = 0.05, h — 0.5, and n — 1, and with the power spectrum P$(fc) normalized to COBE p8| . Although this model 
is almost excluded by current observations, it is still useful to depict the basic effects of the transfer function on the 
bispectrum. 

Figure [l] shows bf{r) (Eq.([l7|)) and bf L {r) (Eq.jl^)) for several different values of r. r = c(tq — r), where r 
is the conformal time, and To is at the present. In our model, ctq = 11.8 Gpc, and the decoupling epoch occurs 
at ct* = 235 Mpc at which the differential visibility has a maximum. Our ctq includes the radiation effect on the 
expansion of universe, otherwise ctq — 12.0 Gpc. r» is the epoch when the most of the primary signal is generated. 
bf (r) and C; look very similar one another in the shape and the amplitude at I > 100, although the amplitude 
in the Sachs- Wolfe regime is different by a factor of —3. This is because Ci is proportional to P$(fc)g|,j(fc), while 
bf{r) cx P$(k)gTi(k), where gn = —1/3. bf (r) has a good phase coherence over wide range of r, while the phase 
of bf L (r) in high-Z regime oscillates rapidly as a function of r. This strongly damps the integrated result of the 
bispectrum (Eq.(|l6|)) in high-Z regime. The main difference between C/ and 6;(r) is that bi(r) changes a sign, while 
Ci does not. 

Looking at figure |l|, we find l 2 b\ ~ 2 x 10 -9 and b^ L ~ 10~ 10 Mpc -3 . The most signal coming from the 
decoupling, the volume element at t* is r 2 Ar* ~ (10 4 ) 2 x 10 2 Mpc 3 , and thus we estimate an order of magnitude of 
the primary reduced bispectrum (Eq.([l9|)) as 



iprimary j— 4 



2r 2 Ar* {l 2 bffbf L x 3 ~T 4 x2x lir 17 f NL . (20) 



Since bf L f N x L ~ r„ 2 <5(r — r„) (see Eq.(p3[)), r 2 /S.r^b^ L f N X L ~ 1. This rough estimate agrees with the numerical result 
below (figure 2). 




FIG. 2. The primary CMB angular bispectrum (Eq.(^)) divided by the Gaunt integral £/™ 
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(Eq.(@)). 
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while the lower panel 



Mi + + 1) (af imi af 2m2 a^ 3 ) /"£ (^ am ») 1 /{^f 

- (9, 11), (99, 101), (199, 201), and (499, 501). 
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Those are shown as functions of I3 for 
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Figure || shows the integrated bispectrum (Eq.(|l6|)) divided by the Gaunt integral !?™/ 2 ™ 2m3 , which is basically 
"hMs '■ Since the signal comes primarily from the decoupling epoch r* as mentioned above, the integration boundary 
is chosen as c (r — 2t*) < r < c (t — O.It*). We use a step-size of O.Ict*, as we found that a step size of O.Olcr* 
gives very similar results. While the bispectrum is a 3-d function, we show different 1-d slices of the bispectrum in 

this figure. h(l 2 + l)h(h + 1) (Kmi a hTn 2 a h™ 3 ) (CiA™ 3 )" 1 / ( 2lr ) 2 is P lotted as a function of ; 3 in the upper panel, 

while ii(ii + l)i 2 (Z 2 + l) K mi «f 2m2 <m 3 ) {QT.IZ T 1 /( 27r ) 2 is P lotted in the lower P aneL i(i + 1 )/( 2 'r) is multiplied 
for each bf (r) which contains P$(fc) so as the Sachs- Wolfe plateau at Z3 < 10 is easily seen in figure 0. l\ and Z2 
are chosen so as (l u l 2 ) = (9, 11), (99, 101), (199, 201), and (499, 501). We find that the (l u l 2 ) = (199, 201) mode, the 
first acoustic peak mode, has the largest signal in this family of parameters. The upper panel has a prominent first 
acoustic peak, and strongly damped oscillations in high-Z regime. The lower panel also has a first peak, but damps 
more slowly. The typical amplitude of the reduced bispectrum is l i bfu %mary f^ L ~ 10~ 17 , which agrees with an order 
of magnitude estimate (Eq.(pfj|)). 

Our formula (Eq.(|nJ)) and numerical results agree with Gangui et al. jl2| calculation in the Sachs- Wolfe regime, 
where gTi{k) ~ —ji(kr*)/3, and thus 

yprirnary ra _ 6/wi (c™ + C?™ Cf 3 W + C™) (Sachs-Wolfe approximation). (21) 

Each term is in the same order as equation (|l^). Cf w is the CMB angular power spectrum in the Sachs- Wolfe 
approximation, 

2 f°° 

Cf w = —J^ fc 2 dfcP $ (fc)jf(fcr,). (22) 
In deriving equation (^) from ©, we approximated bf L {r) (Eq.©) to 

b? L (r) « (-^f 1 ) \ £ *?dkMkr.)Mkr) = -^r; 2 <5(r - r„). (23) 

The Sachs- Wolfe approximation (Eq.(pl|)) is valid only when l\, I2, and I3 are all less than ~ 10, where Gangui et 
al.'s formula gives ~ —6 x 10~ 20 in figure g. It should be stressed again that the Sachs- Wolfe approximation gives the 
qualitatively different result from our full calculation (Eq.(|l9|)) at U > 10. The full bispectrum does change a sign, 
while the approximation never changes a sign because of the use of C ; . The acoustic oscillation and the sign change 
are actually great advantages, when we try to separate the primary bispectrum from various secondary bispectra. We 
shall study this point later. 



C. Primary skewness 

The skewness S3, 



is the simplest statistic characterizing the non-Gaussianity. S3 is expanded in terms of Bi ± i 2 i 3 (Eq.(|3[)) or bi t 
(Eq.@)as 



<? - 1 V J ffii + l) (2k + l) (2Z 3 + 1) ( h h h \ o www 

S " - 4^ V I J B hhh w hW l2 w l3 

= i E + (fa + 5) (is + 5) (0 q) *hh*W h W h W h , (25) 

where Wi is the experimental window function. We have used equation @ to replace Bi^i 2 i 3 by the reduced bispectrum 
bi-tizis in the last equality. Since I = and 1 modes are not observable, we have excluded them from the summation. 

Throughout this paper, we consider the single-beam window function, Wi = e~ l ( l+1 W( 2a ' b \ where Ob = FWHM/V81n2. 
Since ( q q q ) ^1(2(3 1S symmetric under permutation of indices, it is useful to change the way of summation as 



G 



E — > 6 E • ( 26 ) 

2<hhl 3 2<l 1 <l 2 <h 

Since this reduces the number of summations by a factor of ~ 6, we shall use this convention henceforth. 
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FIG. 3. The upper panel shows the primary CMB skewness (Eq.(^)) summed up to a certain 1$, —Ss{< h)f^ x 10 15 . The 
lower panel shows the noise (Eq.(|6o|)) summed up to £3, crs 3 (< I3) x 10 15 . Solid line represents the zero-noise ideal experiment, 
while dotted lines show COBE, MAP, and Planck experiments. 

The upper panel of figure ^| plots Ss(< which is S3 summed up to a certain I3, for FWHM beam-sizes of 7°, 13', 
and 5'.5. These values correspond to beam-sizes of COBE, MAP, and Planck experiments, respectively. Figure || also 
plots the infinitesimal thin-beam case. MAP, Planck, and the ideal experiments measure very similar S3 one another, 
despite the fact that Planck and the ideal experiments can use much more number of modes than MAP. The reason 
is as follows. Looking at equation (p5[), one finds that S3 is the linear integration of b^i^ over i,. Thus, integrating 
oscillations in tf"^ Lrv around zero (see figure ||) damps the non-Gaussian signal in small angular scales, I > 300. 
Since the COBE scale is basically dominated by the Sachs- Wolfe effect, no oscillation, the cancellation affects S3 less 
significantly than in MAP and Planck scales, while Planck suffers from severe cancellation in small angular scales. 
Even Planck and the ideal experiments measure the same amount of S3 as MAP does. As a result, measured S3 
almost saturates at the MAP resolution scale, I ~ 500. 

We conclude this section by noting that when we can calculate the expected form of the bispectrum, then it is a 
"matched filter" for detecting the non-Gaussianity in the data, and thus much more powerful tool than the skewness 
in which the information is lost through the coarse-graining. 

IV. SECONDARY SOURCES OF THE CMB BISPECTRUM 

Even if the CMB bispectrum were significantly detected in the CMB map, the origin would not necessarily be 
primordial, but rather would be various foregrounds such as the Sunyaev-Zel'dovich effect |^9| (hereafter SZ), the 
weak lensing effect, extragalactic radio sources, and so on. In order to isolate the primordial origin from others, we 
have to know the accurate form of bispectra produced by the foregrounds. 
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A. Coupling between the weak lensing and the Sunyaev— Zel'dovich effects 



The coupling between the SZ and the weak lensing effects would produce an observable effect in the bispectrum 
l|,p2[. The CMB temperature field including the SZ and the lensing effects is expanded as 

AT(n) _ AT P (n + V6(n)) , AT sz (n) ^ AT p (n) , y / AT p (n) \ &) + AT 5Z (n) ^ 



/yi I I V / r^i 

where P denotes the primary anisotropy, 0(n) is the lensing potential: 

e(n) = -2/ dr^— ^$(r,nr), (28) 
Jo rr* 



and S*Z denotes the SZ effect: 

AT sz (n) 



y{n)j v , (29) 



T 

where j v is the spectral function of the SZ effect |2{|. y(n) is the Compton y-parameter given by 

drT p (r,hr) _ 2 



y(h)^y Q I - »IJ ' a ~», (.10) 

where 



VTpgartkBTpon _ ln -4, -1 /o »,2\ AijTpO 



^ ^ = 4.3 x 10 -4 -i n . (31) 

p e m p m e c z x V 1 kev / \1" Gpc J 

T p = p gas T e /~p gag is the electron temperature weighted by the gas mass density, the overline denotes the volume 
average, and the subscript means the present epoch. We adopt fi' 1 = 0.88, where \i~ = n e /(p gas /m p ) is the 
number of electrons per proton mass in fully ionized medium. Other quantities have their usual meanings. 
Transforming equation ( p7| ) into harmonic space, 

P , / -] \rn^-mm' m " [0* + l ( l + jj + + jj „P o , „SZ 

a; m - a im + 2^ 2^ V- 1 ) ^ii'i" 2 H'm'Vv'm" + a lm 

I'm' I" m" 



EST ( 1 \m+m'+m" n -mm'm" + + 1) + j(t + 1) Q * / qo \ 

/ 4 yil'l" 9 a J'-m' w l"-m" + a Zm j l^J 

Z'm' I" m" 



where G^l™ 2 " 13 is the Gaunt integral (Eq.(g)). Substituting equation (B2) into (0), and using the identity 



I -I- + I s ! „ . „_ . 1 

(33) 



flgl + 1) - h{h + 1) + / 3 fe + 1) ^ P /Q . SZ \, K 
2 C h \®l 3 m 3 a hm 3 ) + 5 Permutations 



irn2 7713 /imim2m3 



The form of equation (^) is confirmed, and then the reduced bispectrum includes terms in the square bracket. 

The cross-correlation power spectrum of the lensing and the SZ effects, (6* m a^f ), appearing in equation ( |33[ ) was 
first derived by Goldberg and Spergel J2l| . They assumed the linear pressure bias model proposed by Persi et al. |30| : 
T p = T p b gas 5, and the mean temperature evolution of T p ~ T p q(1 + z)^ 1 for z < 2 as roughly suggested by recent 
hydrodynamic simulations 1 31 - 33 1 . Then they derived 

/n* sz\ _ • -iypbgasl 2 f z * dr 2 2 n - r(g) D /, / \ , q/n 



m m H$ J dz VA ' r 2 r5(z) "V K^) / ' 

where £)(2;) is the linear growth factor. Simulations without non-gravitational heating [p2|]33]| suggest that T p o ~ 
0.2 — 0.4 keV and b gas ~ 5 — 10, and similar numbers are obtained by analytic estimations [p2lj34|. In this pressure 
bias model, free parameters except cosmological parameters are T p q and b gas . However, both actually depend on 
cosmological models |§. Since ^ ( @ tm a fm) ~ 2 x 10~ 10 juTpob ga8 @|| and / 2 C ; P - 6 x 10~ 10 , 
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bur lenS ~ ^ 3 W 2 CF) (I 3 (QLAn)) * 5/2] ~r 3 x3x 10" 19 j„T p0 b gas , (35) 

where T p0 is in units of 1 keV, and b hhh = B™^ 2 ™ 3 (G™i™* m3 )~ l (Eq.(§)) is the reduced bispectrum. Thus, 
comparing this to equation (pp[), we obtain 

i primary / p \ 

f& f- 1 x 10 ( 1 NL ). (36) 

lsz — lens \ A T h I 

This estimate suggests that the primary bispectrum is overwhelmed by the SZ-lensing bispectrum in small angular 
scales. This is why we have to separate the primary from the SZ-lensing effect. 



B. Extragalactic radio and infrared sources 

The bispectrum from extragalactic radio and infrared sources whose fluxes F are less than a certain detection 
threshold F d is relatively simple to estimate, when they are assumed to be Poisson distributed. Since the Poisson 
distribution has the white noise spectrum, the reduced bispectrum (Eq.(|^)) is constant, &f x * a | 3 = b ps = constant, then 
we obtain 

where 

F d ) ee g*{x) £ d dFF^ = g^V^JApl (38) 

The assumption of the Poisson distribution is fairly good approximation as found by Toffolatti et al. [p5| . dn/dF is the 
differential source count per unit solid angle, and n(> F d ) = dF(dn/dF). The power law count, dn/dF cx 
with @ < 2, has been assumed, x = hvjk B T ~ (v/56.80 GHz)(T/2.726 K) _1 , and 

= 2 JM1_ /smlW2\ 2 ^ 1 f T \ 3 / sinhx/2 \ 2 

3W " (fc B T)3 ^ a; 2 " 67.55 MJysr-i ^2.726 Ky< \ J ' 1 J 

6 ps is otherwise written in terms of the Poisson angular power spectrum C ps : 

C{< F d ) ee , 9 2 (*) ^ ^F 2 ^ = 5 2 (-)4^^^' ( 4 °) 

as 

V»« F d ) = (2 3 "_^ /2 [n(> F d )]- 1/2 [C^(< F d )] 3/2 . (41) 

Toffolatti et al. @ estimated n(> F d ) - 300 sr" 1 for F d - 0.2 Jy at 217 GHz. This F d corresponds to 5a 
detection threshold for Planck experiment at 217 GHz. Refregier, Spergel and Herbig |3(J extrapolated Toffolatti et 
al.'s estimation to 94 GHz, and obtained n(> F d ) ~ 7 sr _1 for Fd ~ 2 Jy, which corresponds to MAP 5a threshold. 
These values yield 

C ps (90 GHz, < 2 Jy) - 2 x 10~ 16 , (42) 
C ps (217 GHz, < 0.2 Jy) ~ 1 x 10~ 17 . (43) 

Thus, rough estimates for b ps are 

6 ps (90 GHz, < 2 Jy) - 2 x 10" 25 , (44) 
6 ps (217 GHz, < 0.2 Jy) - 5 x 10" 28 . (45) 

While we assumed the Euclidean source count ((3 — 3/2) here for defmiteness, this does not affect an order of magnitude 
estimates above. Since the primary reduced bispectrum cx l~ A (Eq. J20|)) and the SZ-lensing reduced bispectrum cx 
(Eq.(§5j)), the Poisson bispectrum rapidly becomes to dominate the total bispectrum in small angular scales, 

7 primary / f \ 

■Lsz — lens / ■ rp i \ 

mi r 3 x io 6 1 3v p0 9as i (47) 

For example, the SZ-lensing bispectrum measured by MAP experiment is overwhelmed by point sources at I > 100. 
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V. MEASURING BISPECTRA 



A. Fisher matrix 



We shall discuss the detectability of CMB experiments to the primary non-Gaussianity in the bispectrum. We 
also need to separate it from secondary bispectra. Suppose that we try to fit the observed bispectrum by 
theoretically calculated bispectra which include both primary and secondary sources. Then we minimize x 2 defined 
by 



E 



j^obs 

1 1 1 2 1 3 



2<h<l 2 <h 



'hhh 



(48) 



where i denotes a component such as the primary, the SZ and lensing effects, extragalactic sources, and so on. 
Unobservable modes I = and 1 are removed. In case that the non-Gaussianity is small, the cosmic variance of the 



bispectrum is given by the six-point function of ai r 



The variance of -B/ 1 ; 2 ; 3 is then calculated as Pd,p9| 



(49) 



where /\i 1 i 2 i 3 takes values 1, 2, and 6 for cases of that all Vs are different, two of them are same, and all are same, 
respectively. C; = C; + Cf 1 is the total CMB angular power spectrum, which includes the power spectrum of the 



detector noise C, 



N 



is calculated analytically using the formula derived by Knox |4fJ with the noise characteristics 



of the relevant experiments. We do not include C/ from secondary sources, as they are totally subdominant compared 
with the primary Ci and for relevant experiments. For example, inclusion of Ci from extragalactic sources 
(Eqs.([l2"|) or ((43|)) changes our results less than 10%. 
Taking d^cJdAi — 0, we obtain the normal equation, 



E 

3 



E 



2<h<l 2 <h 

Thus, we define the Fisher matrix Fij as 



^l\l 2 lz 



A, = 



V ij — E 

2<ll<l 2 <h 



o(i) dO') 



'hhh 



7T 



2<h<l 2 <h 



E 

2<h<l 2 <l 3 



lihh hhh 



B? b f, B 



hhh 



h + =• )\h + - is + t: 



h h h 




2 h H) b U) 

hhl-3 hhl-3 



hhl-3 



(50) 



(51) 



where we have used equation (^|) to replace Bi ± i 2 i 3 by the reduced bispectrum bi 1 i 2 i 3 (see Eq.(|^) for definition). Since 
the covariance matrix of Ai is i^J , we define the signal-to-noise ratio (S/N)i for a component i, the correlation 
coefficient r,j between different components i and j, and the degradation parameter di of (S/N)i due to as 



F: 



F^FZ 1 



A. = TT.p- 



(52) 



(53) 



(54) 



Note that m does not depend on amplitudes of bispectra, but shapes, di is defined so as di = 1 for zero degradation, 
while di > 1 for degraded (S/N)i. Spergel and Goldberg |20| and Cooray and Hu [E2| considered the diagonal 



component of F i - 1 , while we study all components in order to discuss the separatability between various bispectra. 

An order of magnitude estimation of S/N as a function of a certain angular resolution I is possible as follows. Since 

, I I I ^ 2 

the number of modes contributing to S/N increases as l 3 ^ 2 and / 3 







a-. 



1 

3tt" 



/3 /2 x p/2 



I I I 





l 3 b 



III 



{l 2 d) 3 / 2 



0.36x1, we estimate (S/N), ~ (Fa) 1 / 2 
l 5 b% x 4 x 10 12 , (55) 



10 



-10 



where we have used l 2 Ci ~ 6 x 10" 

Table H and || tabulate all components of i 7 ^ and F,^ 1 , respectively. Table III summarizes (S/N)i 
tabulates di in the diagonal, and in the off-diagonal parts. 



while table IV 



B. Measuring primary bispectrum 

Figure shows the numerical results of differential S/N for the primary bispectrum at In £3 interval, 
1/2 

[d(S/N) 2 /dhxh] Jnl, m the upper panel, and (S/N)(< h)fxL, which is S/N summed up to a certain I3, in 
the lower panel. The detector noises have been computed for COBE 4-yr map fill , for MAP 90 GHz channel, 
and for Planck 217 GHz channel, but the effect of limited sky coverage is neglected. Figure || also shows results for 

the ideal experiment with no noise: Cf = 0. Both [d(S/N) 2 /d\nl 3 ] 1/2 and (S/N)(< h) are monotonically increasing 
function with l 3 as roughly oc l 3 up to l 3 ~ 2000 for the ideal experiment. 

Beyond I3 ~ 2000, an enhancement of the damping tail in C; because of the weak lensing effect p2| stops 

[d(S'/iV) 2 /dln;3] 1/2 and then (S/N)(< l 3 ) increasing. This leads to an important constraint on the observation; 

even for the ideal noise- free and the infinitesimal thin-beam experiment, there is an upper limit on the value of 

1/2 

S/N < 0.3 fNL- For a given realistic experiment, [d(S/N) 2 /d\nl 3 ~\ has a maximum at a scale near the beam-size. 



10° 

10- 1 

10" 2 

10- 3 

10° 

10- 1 

10- 2 

10- 3 




ideal (zero — noise) 

realistic 

(S/N)(<L)f" 1 




1000 



FIG. 4. The predictions of the signal-to-noise ratio, S/N, for COBE, MAP, and Planck experiments (see Eq.([32|)). The 
differential S/N at In £3 interval is shown in the upper panel, while the cumulative S/N up to a certain I3 is shown in the 
lower panel. Both are in units of Jnl- Solid line represents the zero-noise ideal experiment, while dotted lines show the 



realistic experiments mentioned above. The total {S/N)f NL 
experiments, respectively. 



are 1.7 x 10" 3 , 5.8 x 10~ 2 , and 0.19 for COBE, MAP, and Planck 



The total (S/N)f„\ are 1.7 x 10~ 3 , 5.8 x 10~ 2 , and 0.19 for COBE, MAP and Planck experiments, respectively (see 
table III). In order to obtain S/N > 1, therefore, we need /jvl > 600,20, and 5 for each corresponding experiment, 
while the ideal experiment requires Jml > 3 (see table |v|). These values are also roughly obtained by substituting 
equation (p0[) into (p5|), 



I x W- 4 f NL . 



(56) 



primary 



The de grad ation parameters d pr imary are 1.46, 1.01, and 1.00 for COBE, MAP, and Planck experiments, respectively 
(see table IV). This means that MAP and Planck experiments will separate the primary bispectrum from others at 1% 
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or better accuracies. Since the primary and other secondary sources change monotonically in the COBE angular scales, 
COBE cannot discriminate between them very well. In the MAP and Planck scales, however, the primary bispectrum 
starts oscillating around zero, and then is well separated in shape from other secondaries, as the secondaries do not 
oscillate. This is good news for the forthcoming high angular resolution CMB experiments. 



C. Measuring secondary bispectra 

The signal-to-noises for measuring the SZ-lcnsing bispectrum (S/N) sz _i ens in units of \j v \ T p0 b gas are 1.8 x 1CT 



0.34, and 6.2 for COBE, MAP, and Planck experiments, respectively (see table III). T p q is in units of 1 keV. Using 



equations (55) and (p5|), we roughly estimate (S/N) sz _i ens as 

S 



N 



I 2 x W' 6 \j v \T p0 b gas . (57) 



Thus, (S/N) sz -i ens increases with the angular resolution more rapidly than the primary bispectrum (see Eq.(p6[)). 
Since \ju\T p ob gas should be an order of unity, COBE and MAP would not be expected to detect the SZ-lensing 
bispectrum; however, Planck would be sensitive enough to detect, depending on the frequency i.e., a value of j„. For 
example, 217 GHz is totally insensitive to the SZ effect as j v ~ 0, while j„ = —2 in the Rayleigh-Jeans regime. 

The degradation parameters d sz _i ens are 3.89, 1.16, and 1.00 for COBE, MAP, and Planck experiments, respectively 
(see table [V). Thus, Planck will separate the SZ-lcnsing bispectrum from other effects. Note that {S/N) sz _i ens 



values must be an order of magnitude estimation, as our cosmological model is the COBE normalized SCDM yielding 
(78 = 1-2. Since this as is about a factor of 2 greater than the cluster normalization with Q m = 1, and 20% greater 
than the normalization with Cl m — 0.3 P3fl . Thus, this factor tends to overestimate (®i m af^) (Eq.Q)) by a factor of 
several. On the other hand, using the linear power spectrum for P$(k) rather than the non-linear power spectrum tends 
to underestimate the effect by a factor of several at I ~ 3000 p2| . However, our main goal is to discriminate between 
shapes of various bispectra, not amplitudes, so that this factor does not affect our conclusion on the degradation 
parameters di. 

For the extragalactic radio and infrared sources, we estimated the signal-to-noises as 5.7 x 10~ 7 (y s /10~ 25 ), 



2.2(bP s /10- 25 ), and 52(&p s /10" 27 ) for COBE, MAP, and Planck experiments, respectively (see table [H]), and the 
degradation parameters d ps are 3.45, 1.14, and 1.00 (see table IV). Since 



from equation (|55|), S/N of the bispectrum from point sources increases very rapidly with the angular resolution. Our 
estimate that MAP will detect the bispectrum from point sources is consistent with the results found by Refregier, 
Spergel and Herbig p&. Although MAP cannot separate the Poisson bispectrum from the SZ-lensing bispectrum 



very well (see in table [V), it would not matter as the SZ-lensing bispectrum would be too small to be measured 
by MAP. Planck will do an excellent job on separating all kinds of bispectra, at least including the primary signal, 
SZ-lensing coupling, and extragalactic point sources, on the basis of the shape difference. 



D. Measuring primary skewness 

For the skewness, we define S/N as 



\N 

where the variance is Hi 



(59) 



2 ' 

a s 3 



-1 

v (2*1 - 1) (2Z 2 + 1) (2l 3 + 1) f h h h \ ., r w i w i w i 
6 (4^)3 [ o ) c hCiA,W h W h W h 



2 



2 



i E ('1 + I) (fa + I) ( *3 + I) ( XWiMwlW?, (60) 
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In the last equality, we have used the symmetry of summed quantity with respect to indices (Eq.([26h), and removed 
unobservable modes I — and 1. Typically as 3 ~ 10 -15 , as as 3 ^JC(O)] 3 ^ 2 ~ 10~ 15 , where C(0) is the temperature 
auto correlation function including noise. The lower panel of figure |3J shows os 3 (< Z3), which is as 3 {< h) summed up 
to a certain 1%, for COBE, MAP, and Planck experiments as well as the ideal experiment. Since C\Wf — C;e~ / ( i+1 ) CT '> + 
u> -1 , where u> _1 determines the white noise power spectrum of the detector noise according to the Knox's formula 
po! , the dominance of second term beyond the experimental angular resolution scale, I ~ cr^ 1 , keeps o"s 3 (< ^3) slightly 
increasing with Z 3 , while Sz(< I3) becomes constant beyond that (see the upper panel of figure |J). As a result, S/N 
starts somewhat decreasing beyond the resolution. We use the maximum S/N for estimating the minimum value of 
f NL needed to detect the primary S3. We find that f NL > 800, 80, 70, and 60 for COBE, MAP, Planck, and the 
ideal experiments, respectively, with all-sky coverage. 

These f^L values are systematically larger than those needed to detect Bi ± i 2 i 3 by a factor of 1.3, 4, 14, and 20, 
respectively (see table |v|). Higher the angular resolution, less sensitive measuring the primary S3 than Bi ± i 2 i 3 . This 
is because the cancellation effect in smaller angular scales due to the oscillation of Bi ± i 2 i 3 damps S3. 



1U" 




1 

(S/N)«l a )fij2 






10- 1 




Bispectrum 




Planck- 
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FIG. 5. The comparison of the signal-to-noise ratio summed up to a certain I3, S/N(< for the detection of the 
bispectrum (upper panel; Eq.(j5^)) and the skewness (lower panel; Eq.(S)_) in units of Jnl for COBE, MAP, and Planck 
experiments (dotted lines), and the ideal experiment (solid line). See tableM for values of Jnl in order to obtain S/N > 1. 

VI. DISCUSSION AND CONCLUSION 

Using the full radiation transfer function , we have computed numerically the primary cosmic microwave background 
bispectrum (Eq. (|l6|)_) and skewness (Eq.(p5[)) down to arcminutes angular scales. The primary bispectrum oscillates 
around zero (figure ||), thus the primary skewness saturates at the MAP angular resolution scale, I ~ 500 (figure ||). 
We have introduced the "reduced" bispectrum defined by equation (^), and confirmed that this quantity is more 
useful to describe the physical property of the bispectrum than the full bispectrum (Eq.(^)). 

Figure p| co mpares the expected signal-to-noise ratio for detecting the primary non-Gaussianity based on the bispec- 
trum (Eq.(|52D) to that based on the skewness (Eq.(|5^)). It shows that the bispectrum is almost an order of magnitude 
more sensitive to the non-Gaussianity than the skewness. We conclude that when we can compute the predicted form 
of the bispectrum, it becomes a "matched filter" for detecting the non-Gaussianity in the data, and thus much more 
powerful tool than the skewness. Table |v| summarizes Jnl required for detecting the primary non-Gaussianity using 
the bispectrum or the skewness with COBE, MAP, Planck, and the ideal experiments. This shows that even the ideal 
experiment needs Jnl > 3 in order to detect the primary bispectrum. 
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We estimated the secondary bispectra from the coupling between the Sunyaev-Zel'dovich (SZ) and the weak lensing 
effects, and from the extragalactic radio and infrared sources. Only Planck will detect the SZ-lensing bispectrum, 
while both MAP and Planck will detect the bispectrum from extragalactic point sources (table 111). 

We also studied how well we can discriminate among the primary, the SZ-lensing coupling, and the extragalactic 
point sources bispectra. We found that MAP and Planck will separate the primary from other secondary sources at 
1% or better accuracies. This conclusion is due to the presence of acoustic oscillation in the primary bispectrum that 
does not appear in the secondary bispectra. The SZ-lensing coupling and the extragalactic sources are well separately 
measured by Planck experiment, although COBE and MAP cannot discriminate between them (table [j"v|) . 

Our arguments about the ability to discriminate among various bispectra were fully based upon the shape difference, 
and thus did not take into account the spectral difference in the frequency space. As pointed out by f^Jfql , the multi- 
band observation is so efficient to discriminate among the primary signal and the other foreground contaminants for 
measuring the CMB anisotropy power spectrum. Their scheme should be effective on the bispectrum as well, and the 
accuracy of the foreground removal will be improved further. Thus, we expect that MAP and Planck will measure 
the primary bispectrum separately from the foregrounds. 

The simplest inflationary scenario usually predicts small Jnl (~ 1CP 2 ) fll]Jl^] , and the second order perturbation 
theory yields Jnl ~ 1 Pq)- Thus, the significant detection of the primary bispectrum or the skewness with any 
experiments means that the simplest inflationary scenario needs to be modified. According to our results, if the 
reported detections 0,|l8| of the bispectrum in the COBE map were the cosmological origin, then MAP and Planck 
would detect the primary bispectrum much more significantly. Although Banday, Zaroubi and Gorski ||l9f pointed 
out the one of those detections |T^] could be accounted for by the experimental systematic effects of COBE, the other 
is claimed to be significant even after removing such the systematics. 

Although we have not discussed so far, the spatial distribution of emissions from interstellar dust is a potential 
source of the microwave non-Gaussianity. Since it is very hard to estimate the bispectrum analytically, the dust map 
compiled by Schlegel, Finkbeiner and Davis |47[ | could be used to estimate the dust bispectrum. For example, we 

— 3/2 

found that the dimensionless skewness parameter defined by ((AT) 3 ) /((AT) 2 ) is as large as 51. We used the 
publicly available HEALPix-formatted E§] 100 /im map which contains 12,582,912 pixels without sky cut. The mean 
intensity in the map was 14.8 MJy sr _1 . Of course, this skewness is largely an overestimate for the CMB measurement 
in reality; we need to cut a fraction of sky which contains the Galactic plane, and then this will greatly reduce the 
non-Gaussianity. Nevertheless, residual non-Gaussianity is still a source of the microwave bispectrum, and has to be 
taken into account. Moreover, the form of the bispectrum measured in the dust map would reflect the physics of 
interstellar dust, which is highly uncertain at present, and thus studying the interstellar dust bispectrum would be 
challenging field. 
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TABLE I. 
respectively. T 


The 

p o is 


\ 1 V 

Fisher matrix Fy (Eq.(|5lJ)). 
in units of 1 keV, b% = V s /10 


i and j denote components listed in 
" 25 , and 6f* = fe^/lO" 27 . 


the first row and the first column, 


COBE 




primary 


SZ-lensing 


point sources 


primary 
SZ-lensing 
point sources 




4.2 x lO" 6 f NL 


-4.0 X 10" 7 }NLjvTp0b a as 

1.3 x HT 7 {j v T p0 b gas ) 2 


-1.0 x 10" 9 f NL b p 2 l 
3.1 x 10" 10 j„T p0 b gaa b% 

i.i x nn 12 {bUY 


MAP 


primary 
SZ-lensing 
point sources 




3.4 x 10" a f NL 


2.6 x 10~ 3 f N Lj»T p0 b gas 
0.14 (j„Tp b gas ) 2 


2.4 x 10" 3 f NL b v 2 l 
0.31 j„T pobgasb^l 
5.6 (6^) 2 


Planck 


primary 
SZ-lensing 
point sources 




3.8 x 10" 2 f NL 


7.2 X 10" 2 fNLjuT p0 b g a s 

39 (jvTpobga*) 2 


1.6 x 10" 2 /iVL&f? 
5.7 j v T pobgasb 7 ^ 
2.7 x 10 3 (6^) 2 


TABLE II. 
respectively. T 


The inverse Fisher matrix F^ 1 . 
p o is in units of 1 keV, 6fg = 6 ps /10 


i and j denote components listed in 
" 25 , and 6f* = b^/lO" 27 . 


the first row and the first column, 


COBE 




primary 


SZ-lensing 


point sources 


primary 
SZ-lensing 
point sources 




3.5 x 10 5 f~l 


1.1 X 10 B UNLjuTp b gas )~ L 

3.1 x 10 7 (j v T p0 b gas )- 2 


1.3 X 10 Y ifNLbH)- 1 
-7.8 x 10 9 {jvTpobgasbH)- 1 

3.1 x 10 12 {bH)- 2 


MAP 


primary 
SZ-lensing 
point sources 




3.0 x 10 2 f N 2 L 


— 6.1 (fNLjvTpob gas ) 1 

8.4 (>r p ob 9 a S )- 2 


0.21 tfNLbH)- 1 
-0.46 (juTpobga^)- 1 

0.21 {bll)- 2 


Planck 


primary 
SZ-lensing 
point sources 




26 f N l 


-4.9 X W~ 2 (fNLjvTpobgas)- 1 
2.6 X 10" 2 [j v T pQbgas)- 2 


-5.7 x 10" b (fNLb^y 1 
-5.4 x 10" 5 (jvTpobgasb^)- 1 
3.7 x 10" 4 (6^)- 2 


TABLE III. The signal-to-noise ratio (S/N)i (Eq.(|5^)). i denotes a component listed 
keV, b p 2 l = b ps /10~ 25 , and = b ps /10" 27 '. 


in the first row. T p o is in units of 1 






primary 


SZ-lensing 


point sources 


COBE 
MAP 
Planck 




1.7 x 10~ 3 f NL 

5.8 x 10" 2 f NL 

0.19 fNL 


1.8 X 10" 4 \ju\T p0 bgas 

0.34 \j v \T p0 b gas 

6.2 \jv\Tpdbgas 


5.7 x 10" 7 b p 2 l 
2.2 bll 
52 bf? 
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TABLE IV. The degradation parameter d t (Eq.(|54|)) and correlation ry (Eq.(|53|)) matrix, i and j denote components listed 
in the first row and the first column, respectively, di for i = j, while m for i 7^ j. 



COBE 


primary 


SZ-lensing 


point sources 


primary 


1.46 


0.33 sgn(>) 


1.6 x to -2 


SZ-lensing 




3.89 


-0.79 sgn(>) 


point sources 






3.45 


MAP 


primary 


1.01 


-0.12 sgn(>) 


2.7 x 10 _:i 


SZ-lensing 




1.16 


-0.35 sgn(>) 


point sources 






1.14 


Planck 


primary 


1.00 


-5.9 x 10~ 2 sgn(>) 


-5.8 x 10~ 4 


SZ-lensing 




1.00 


-1.8 x 10~ 2 sgn(>) 


point sources 






1.00 


TABLE V. The minimum non- linear couplin 


g constant }nl required to detect the primary 


non-Gaussianity by the bispec- 


trum and the skewness with the 


signal-to-noise 


ratio of > 1. These estimates include the effects of cosmic variance, detector 


noise, and foreground sources. 








Experiments 




}nl (Bispectrum) 


fNL (Skewness) 


COBE 




600 


800 


MAP 




20 


80 


Planck 




5 


70 


Ideal 




3 


60 
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